We claim: 



CLAIMS 



1 . A method for editing a geometrical model with a level set modeling surface 
editor operator, comprising: 

performing a level set surface editing operation on a level set model, wherein said 
operation is defined by a level set surface editing operator. 

2. The method of claim 1 further comprises: 

converting an input model into said level set model for said step of performing a 
level set surface editing operation; 

rendering said level set model after said step of performing a level set surface editing 
operation. 

3. The method of claim 2 wherein said converting step converts said input 
model from a geometric to a volumetric model by using scan conversion. 

4. The method of claim 2 wherein said converting step converts said input 
model with by using distance calculations. 

5. The method of claim 2 wherein said converting step converts said input 
model with by using Sethian's Fast Marching Method. 
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6. The method of claim 1 wherein said level set model is represented in a 
narrow-band distance volume. 

7. The method of claim 1 wherein said performing a level set surface editing 
5 operator is defined as a speed function, said function comprising: 

a regional constraint function component; 

a filter function component; 

a surface properties defining function component. 

10 8. The method of claim 7 wherein said speed function is 

F(x,n,</>) = D q {d)C(rt)G(rt) , wherein 

D q (d) is said regional constraint function component, 

C(7) is said filter function component, and 

G(7) is said surface properties defining function component. 

15 

9. The method of claim 8 wherein said D q (d) is further defined as having q be 

a geometric primitive such as a point, a line-segment, or surface and d be the distance from 
the level set surface of said level set surface model to q . 

20 10. The method of claim 8 wherein said 7 is a geometric property of the level 

set surface of said level set surface model, wherein 7 is any order property of (f> . 
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11. The method of claim 10 wherein 7 represents distance, which is the zero-th 
order property of <f> . 

12. The method of claim 10 wherein 7 represents normal vector to said level set 
5 surface, which is the first order property of 0 . 

13. The method of claim 10 wherein 7 represents curvature of said level set 
surface, which is the second order property of <f> . 

10 14. The method of claim 8 wherein said regional constraint component defines a 

region-of-influence. 

15. The method of claim 14 wherein said region-of-influence is defined by a 
distance calculation to a geometric primitive such as a point set, a line segment, or a surface. 

15 

16. The method of claim 15 wherein said region-of-influence is defined by a 
distance to an intersection curve point set. 

17. The method of claim 15 wherein said region-of-influence is a super-ellipsoid. 

20 

18. The method of claim 8 wherein said filter function component is a filter of 
local geometric surface properties, wherein user can control the behavior of said level set 
surface editing operator based on a function of a local geometric surface property. 
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19. The method of claim 8 wherein said surface properties defining function 
component defines the behavior of said level set surface editing operator. 

20. The method of claim 1 wherein said level surface editing operator is a CSG 
5 intersection operator. 

21. The method of claim 1 wherein said level surface editing operator is a CSG 
difference operator. 

10 22. The method of claim 1 wherein said level surface editing operator is a CSG 

union operator. 

23. The method of claim 7 wherein said level surface editing operator is a 
blending operator, wherein blending is constrained by a region-of-influence. 

15 

24. The method of claim 23 wherein said region-of-influence is based on the 
distance to an intersection curve shared by both input surfaces. 

25. The method of claim 23 wherein said intersection curve is represented by a 
20 point set. 

26. The method of claim 23 wherein said blending operator is defined by the 
function F bhnd (x, n, (f>) = aD q (d)C(K)K , wherein 
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a is a user-defined positive scalar that controls the rate of the level set calculation; 
D q (d) is said region-of-influence component with d being the shortest distance 

from the level set surface to said intersection curve point set; 
K is a curvature term; and 

C(K) is a filtering function that defines the geometric properties of said blending 
operator. 

27. The method of claim 23 wherein said blending operator calculates a closest 
point in a set by using the approximate nearest neighbor search algorithm of Mount and 
Arya. 

28. The method of claim 27 wherein said blending operator calculates a closest 
point in a set by storing a point set in a K-D tree, wherein points in said point set is are 
uniform distributed. 

29. The method of claim 27 wherein said blending operator calculates a closest 
point in a set by storing a point set in the balanced box decomposition (BBD) tree, wherein 
points in said point set is are clustered. 

30. The method of claim 7 wherein said level surface editing operator is a 
smoothing operator, wherein a surface is smoothed by applying motions in a direction that 
reduces local curvature. 
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31. The method of claim 30 wherein said smoothing operator is constrained to 
move outward relative to said surface to smooth said surface by adding material to said 
surface. 

32. The method of claim 30 wherein said smoothing operator is constrained to 
move inward relative to said surface to smooth said surface by removing material to said 
surface. 

33. The method of claim 30 wherein said smoothing operator is converted into a 
sharpening operator by changing the sign of a scalar term and using said filter function to 
truncate said speed-function to zero for high values of the curvature. 

34. The method of claim 30 wherein said smoothing operator is defined by the 
function F smooth (x,n 9 <f>) = aD s (d)C(K)K , wherein 

a is a user-defined positive scalar that controls the rate of the level set calculation; 

D s (d) ensures that said function smoothly goes to zero near the boundary of the 

region-of-influence; 

A" is a curvature term; and 

C(K) is a filtering function that defines the geometric properties of said smoothing 
operator. 
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35. The method of claim 34 wherein said region-of-in£luence of said D s (d) is 

defined said by a distance calculation to a geometric primitive such as a point set, a line 
segment, or a surface. 

36. The method of claim 35 wherein said smoothing operator calculates a closest 
point in a set by using the approximate nearest neighbor search algorithm of Mount and 
Arya. 

37. The method of claim 36 wherein said smoothing operator calculates a closest 
point in a set by storing a point set in a K-D tree, wherein points in said point set is are 
uniform distributed. 

38. The method of claim 36 wherein said smoothing operator calculates a closest 
point in a set by storing a point set in the balanced box decomposition (BBD) tree, wherein 
points in said point set is are clustered. 

39. [new] The method of claim 35 wherein said region-of-influence is a super- 
ellipsoid. 

40. The method of claim 7 wherein said level surface editing operator is a point 
set attraction/repulsion operator. 
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41. The method of claim 40 wherein said point set attraction/repulsion operator 
calculates a closest point in a set by using the approximate nearest neighbor search algorithm 
of Mount and Arya 

42. The method of claim 41 wherein said point set attraction/repulsion operator 
calculates a closest point in a set by storing a point set in a K-D tree, wherein points in said 
point set is are uniform distributed. 

43. The method of claim 41 wherein said point set attraction/repulsion operator 
calculates a closest point in a set by storing a point set in the balanced box decomposition 
(BBD) tree, wherein points in said point set is are clustered. 

44. The method of claim 40 wherein said point set attraction/ repulsion operator 
is used to emboss a surface. 

45. The method of claim 1 wherein said level surface editing operator is a 
morphological editing operator. 

46. The method of claim 2 wherein said rendering uses a volume rendering. 

47. The method of claim 2 wherein said rendering uses Incremental Marching 
Cubes mesh extraction, which optimizes the mesh extraction by limiting the region which 
needs to be processed. 
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48. The method of claim 2 wherein said input model is a polygon mesh. 



49. The method of claim 48 wherein said polygon mesh is scan converted into a 
level set model by computing a distance volume, wherein said computing calculates the 
closest point on and shortest singed distance to said mesh by solving the Eikonal equation 

|V0| = 1. 

50. The method of claim 49 wherein said computing a distance volume uses the 
CPT (closest point) algorithm of Mauch. 

51. The method of claim 2 wherein said input model is a Constructive Solid 
Geometry (CSG) model. 

52. The method of claim 2 wherein said input model is an implicit model. 

53. The method of claim 2 wherein said input model is a scanned volume. 

54. The method of claim 8 further comprises: 

resetting the volumetric representation of said level model after said step of 
performing a level set surface editing operation to ensure that <f> is approximately equal to 
the shortest distance to the zero level set in the narrow band. 
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55. The method of claim 7 wherein level set computations associated with 
solving said speed function are regionally constrained to a sub-volume defined by the 
bounding-box of a region-of-influence primitive. 
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